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Abstract. We analyze the possibility that the spatial pe- 
riodicity of 128/i _1 Mpc in the galaxy count number ob- 
served recently in deep pencil-beam surveys could be ex- 
plained by the mere existence of structure with the ap- 
propriate scale in the galaxy distribution in the Universe. 
We simulate a universe where the distribution of galax- 
ies has an intrinsic length scale and then investigate the 
probability of observing in a given direction a spatial pe- 
riodicity similar to that observed in the pencil-beam sur- 

: veys of the galactic polar regions. A statistical analysis 
shows that this probability is of the order of 10 -8 , a value 
which excludes the afore mentioned explanation. This re- 
sult contrasts with the estimates based on Voronoi foam 

' models which yield probabilities of the order 10 _1 — 10 -2 . 
The conclusion that emerges from these studies put to- 
gether with the present one, is that although they can 

, be obtained in some models, the presence of an intrinsic 
scale in the structure is not sufficient to achieve reasonable 
probabilities for observation of the detected periodicities. 
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1. Introduction 

The cosmic microwave background radiation shows a re- 
markable isotropic structure of the order ST/T ss 10~ 5 , 
which is taken as an evidence for the smoothness of the 
mass density of the Universe on the largest scales com- 
parable to the Hubble radius (Wilkinson 1987), Hq 1 w 
3000ft- 1 Mpc {H Q = 100ft km s" 1 Mpc" 1 ). On the other 
hand, the existence of galaxies and clusters of galaxies as 
well as the (larger than unity) amplitude of the galaxy- 
galaxy correlation function (Peebles 1980) indicate that 
the Universe possesses a highly irregular structure on 
scales less than 10ft. -1 Mpc. 
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On intermediate scales, ranging from about 30ft- 1 Mpc 
to around 2000ft -1 Mpc, the existing data do not allow a 
conclusive interpretation about the mass density distri- 
bution of the Universe. Nevertheless, there are a num- 
ber of different observations suggesting that the structure 
on these scales is quite varied. The most interesting of 
these observations include the voids with sizes greater 
than 30ft -1 Mpc observed in the distribution of bright 
galaxies (Kirschner et al. 1981; de Lapparent et al. 1986) 
the (larger than unity) amplitude of the cluster-cluster 
correlation function on scales up to 30ft" 1 Mpc, the pecu- 
liar velocity field in our local vicinity on a scale of about 
50ft -1 Mpc, and the "Great Wall" at a distance of ap- 
proximately 100ft -1 Mpc (Geller & Huchra 1980). Also, 
noteworthy is the analysis of the spatial distribution of 
superclusters on the intermediate scales (Einasto et al. 
1997). For this analysis a sample of more than 200 super- 
clusters (systems with more than 8 clusters where the dis- 
tances between nearest-neighbors among member clusters 
do not exceed 24ft -1 Mpc) was used, where the estimated 
redshift of simple members reaches up to z = 0.12 (around 
1 /4 of the redshift reached in the pencil-beam galactic sur- 
veys). These three-dimensional data of the distribution of 
superclusters has been used to calculate the correlation 
function as well as the power spectrum. In the analysis of 
the power spectrum it was found that it is characterized 
by a peak with a wavelength of (120 ± 15)ft -1 Mpc. 

But probably the most puzzling discovery concerns the 
recent observations which show an apparent periodicity 
in the distribution of redshifts in pencil-beam surveys of 
the north and south galactic polar regions (Broadhurst et 
al. 1990). These surveys extend to a redshift of around 
z w 0.5 in both directions, and there is a strong evidence 
for a periodicity of 128ft -1 Mpc which extends for over 13 
periods. Although this result was considered initially as 
a statistical anomaly (Kaiser & Peacock 1991), later and 
more detailed studies (Kopylov et al. 1988; Mo et al. 1992; 
Fetisova et al. 1993; Einasto & Gramann 1993; Szalay et 
al. 1993; Einasto et al. 1997) provided new evidence for 
the alluded periodicity on the same scales. 
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If this periodicity would be established in most di- 
rections, a naive interpretation would imply that galax- 
ies and clusters are distributed on concentric shells cen- 
tered around our galaxy, clearly a blow to our cosmological 
conceptions. Another more plausible explanation is based 
upon the idea that the observed periodicity is just a vi- 
sual effect induced by a periodic oscillation of the gravita- 
tional constant (Morikawa 1990; Hill, Steinhardt & Turner 
1990; Salgado, Quevedo & Sudarsky 1996, 1997; Quevedo, 
Salgado & Sudarsky 1997; Sudarsky, Quevedo & Salgado 
1998), and that the mass density of the Universe on inter- 
mediate scales remains homogeneous and isotropic as on 
the largest scales. This explanation requires the introduc- 
tion of a time-dependent scalar field into the model (the 
oscillating G model) which is non-minimally coupled to 
gravity. 

In view of the remarkable conclusions that would de- 
rive from any of these explanations it seems imperative 
to consider whether the observations are likely to occur 
in the absence of these scenarios, and can be the result of 
chance augmented with the assumption of the presence of 
structure at the appropriate scale but with an otherwise 
random distribution. We will study the plausibility of this 
latter explanation inspired on the idea, suggested by Sza- 
lay et al. (1993), that the structure of the Universe on 
these scales consists of a regular network of superclusters 
and voids with a step size of the order of the wavelength 
discovered in the power spectrum. The specific structure 
that is studied by Van de Weygaert (1991) and Subba Rao 
and Szalay (1992) resembles a collection of honeycombs or 
a three-dimensional chess-board. Such structures known 
as Voronoi foams result in a reasonable probability for the 
observation of a periodicity similar to that in the actual 
observations. 

The aim of this work is to investigate whether it is 
sufficient to have a universe with an intrinsic length scale 
in its large scale structure in order to explain the observa- 
tion of the periodicity mentioned above. In particular, the 
present work can be seen as a test of the suggestion that a 
hypothetical "built-in scale in the initial spectrum" arising 
from double inflation models and supported by cosmic mi- 
crowave background radiation data (Einasto, 1997), could 
be the explanation of the above mentioned periodicity 
Therefore, we will consider a model that does not intro- 
duce any other large scale pattern besides the existence of 
the above mentioned characteristic scale. 

Needless is to say that neither our model nor the 
Voronoi foam models are anything more than "toy" mod- 
els as they lack, for example, a proven dynamical scenario 
for their formation starting from reasonable initial data. 
We want to construct a simple realization of such a model. 
We will do this by assuming that the galaxies are evenly 
distributed over the surface of three-dimensional spheres 
with an average radius comparable with the appropriate 
scale, with the spheres themselves randomly distributed 
in the Universe. In a universe filled by these spheres we 
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Fig. 1. Two-dimensional section of a universe filled with 
a random smearing of three-spheres with a Gaussian dis- 
tribution in radii. A constraint is imposed in order that 
the spheres do not intersect with each other. Notice that 
the distortion of the circles into ellipses is only due to the 
difference of scales on the x— y axes. 

choose a point as the "observational site" and perform 
a statistical analysis to investigate the probability of ob- 
serving a periodic distribution with the characteristics of 
those reported in (Broadhurst et al. 1990) in a particular 
direction. We obtain an upper bound for the value of such 
probability of the order 10~ 8 , a value that obviously ex- 
cludes the alluded structure as an acceptable explanation 
of the observations. The assumption that the galaxies are 
homogeneously distributed on the surface of the spheres 
is not expected to affect the main result of our numerical 
analysis in the direction that would invalidate this upper 
bound. 

2. The model 

We have taken the radii of the spheres to have a Gaus- 
sian distribution centered around the mean value of 80/i _1 
Mpc and analyzed the results for different choices of the 
standard deviation a of the distribution (namely with 
(j = 5, 15,30/i _1 Mpc). We have implemented a code 
that simulates the universe as a three-dimensional vol- 
ume of (5000 h~ 1 Mpc) 3 filled by these spheres subjected 
to the condition that the spheres do not intersect with 
each other. In figure 1, a two-dimensional section of this 
volume is depicted. 

We emphasize that the above condition of no-sphere 
overlapping and the Gaussian distribution on the spheres 
radii imply that the numerical process of filling the toy 
universe stalls after a certain number of spheres. After 
that stage, the inclusion of more spheres demands a very 
long computing time. The adequate number of spheres 
needed to fill the universe is not fixed in advance, but the 
important point is that such a number must be so that the 
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correlation function of the matter in the universe exhibits 
a characteristic scale, which is the most important aspect 
in the scenario that we are testing. We have confirmed 
that once correlations arc present the results become in- 
sensitive to a further increase in the number of spheres. 

In order to check the appearance of scale, we have com- 
puted a two-point correlation function defined in the fol- 
lowing way. Let N be the total number of spheres in the 
toy universe, and let rij(r) be the number of sphere centers 
contained between two concentric shells with radii r and 
r + Sr around the center i (with Sr <~ 10 Mpc h~ l ). The 
number density of sphere centers within r and r + Sr is 
then given by 

This is to be constant if correlations between neighbors 
at all scales are absent. It is customary to extract from 
this, the excess density number by defining the two-point 
correlation function as the difference between the above 
density distribution and the mean density of sphere-center 
number: 
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Fig. 2. Two-point correlation function of the matter con- 
tent as a function of the distance. The dotted, solid and 
dash-dotted lines corresponds to the values a = 5, a = 30 
and a = 15 respectively. Notice the presence of correlation 
at scales of ~160-180 Mpc and anticorrelations at shorter 
scales. 



£(r) = C(r)- < C > . 



(2) 



More specifically, < C > corresponds to the mean of C (r) 
at scales for which C(r) is nearly constant (i.e., £ ~ at 
those scales). If the Universe has a characteristic length 
scale we expect it will be evident in the plot of the corre- 
lation function. 

The results are shown in figure 2 with a typical num- 
ber TV w 4 x 10 4 . We can appreciate the presence of a 
characteristic scale at distances which are comparable to 
the mean value of the spheres radii, and stress that such 
scale is already present when the spheres fill a volume of 
only about 25% of the total volume of the toy universe. 
We do not expect any drastic changes in the qualitative 
as well as in the quantitative behavior of the correlation 
function if we increase the number of spheres in the uni- 
verse. Moreover, as shown in figure 2, the resulting length 
scale is almost independent of the value of a considered. 

Now we will estimate the probability of observing in a 
specific direction a periodicity like the one seen in the real 
observations. To this end we note that the power spectrum 
of the observations in Szalay ct al. (1993) revealed a peak 
at a frequency corresponding to 128/i _1 Mpc which con- 
tained 23% of the total integrated power. Thus we need to 
estimate the probability of observing in a specific direction 
a matter distribution with a power spectrum containing a 
peak with approximately the same percentage of the total 
power. 

The "observational site" has been chosen to be a vertex 
of the box which in fact represents an octant of a larger 
toy-universe. This particular choice allows to consider the 



longest possible extent for the line of sight in each di- 
rection, but it does not affect the result of the statistical 
analysis which is performed for about 10 5 directions and 
for "pencil-beams" extending approximately to 5000/i _1 
Mpc. 

In this model each direction corresponds to a numer- 
ical pencil-beam, that is, a straight line which intersects 
several spheres, each intersection thus representing high- 
density spots. The ensemble of the distance separations 
between the spots results thus in a one-dimensional dis- 
tribution of matter along the line of sight, which we repre- 
sented mathematically by an ensemble of step functions of 
width 10 Mpc h^ 1 with height normalized to unity cen- 
tered on each intersection of the line of sight with the 
surface of a sphere. We have also consider the effect of re- 
placing step functions by Gaussian functions and showed 
that this does not change the probability estimates be- 
yond the statistical error bars. Finally, we have performed 
a Fourier analysis of this distribution and confronted the 
resulting power spectrum with the observational one. Fig- 
ure 3 shows the resulting power spectrum of the matter 
distribution in a representative direction. We note that 
there is a substantial power at high frequencies, a feature 
associated with the use of step functions that is slightly 
suppressed by the use of Gaussian functions. 

In order to compare with the real observations we 
look in each direction for the peak containing most of the 
power. We then calculate the percentage II of the power- 
content within the peak: 



„ Power content withm the peak 

n = — x ioo . 

Total power 



(3) 
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Table 1. Parameters of the power-law (Q) for the proba- 
bility of finding a spatial periodicity in a given direction. 
Here the statistical uncertainties of data are taken into 
account for determining the ranges of the parameters. 
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Fig. 3. Example of a power spectrum (Fourier analysis) of 
the resulting distribution of matter along a line of sight. 



o 




50 100 150 200 



X[Mpc] 

Fig. 4. Distribution of the number of peaks with a given II 
(maximum relative power content) associated to each di- 
rection of sight as a function of their corresponding wave- 
length in the Fourier space. The dotted, solid and dash- 
dotted lines corresponds to the values a = 5, a — 30 and 
a = 15 respectively. 

As we mentioned, this procedure has been applied to 10 5 
different directions starting from the specific site described 
above, obtaining representative peaks with values II at 
different frequencies in each case. The distribution of the 
number of peaks with a given II is plotted in figure 4 
as a function of its corresponding wavelength. We note 
that most of the highest peaks cluster around a specific 
wavelength of ~ 40 Mpc h _1 , however the power contained 
in most of them represents barely 2% of the corresponding 
total integrated power (see for instance the figure 3). 

Since one can define the probability of observing a pe- 
riodicity in a certain direction as the number of directions 
in which the periodicity is present over the total num- 



a c ax 2 

30 2.446564 x 10 6 ± 1.722128 x 10 6 -7.30 ± 0.20 0.401 

15 2.390874 x 10" ± 1.682928 x 10 b -7.22 ± 0.20 0.377 

5 1.3196748 x 10 Y ± 1.0245539 x 10' -8.81 ± 0.31 2.468 



ber of directions, one should analyze an adequate sample 
in which one observes a sufficient number of cases with 
the desired periodicity. Since the sample of 10 5 directions 
showed no periodicity with a peak having the required 
power II ~ 23%, one should be tempted to give only an 
upper bound of 10~ J for the probability of observing that 
value. This is not the best upper bound that can be ob- 
tained and it is only determined by the finiteness of the 
computing time. We have been able nevertheless to im- 
prove on this bound by the following procedure: We have 
fitted our data for the probability p of a peak with a given 
percentage II of the power as a function of that percent- 
age to various analytical expressions and found that this 
function is very well fitted by a power law: 

p(n) = cu a (4) 

where c and a are constants with values depending on 
the chosen a (see table 1). This allows us to extrapolate 
the probability of finding the observed periodicity with 
II = 23% in a given direction. 

An example of this law for a = 30 Mpc is shown in 
figure 5. A similar fit was made for the cases with other 
values of a with the remarkable result that the power law 
is practically unaffected. This is a strong evidence for the 
robustness of our results. 

We also stress that the data-set fitted to the above 
power law has values of x 2 , which according to the degrees 
of freedom of the data (~ 2), indicates the good quality 
of the fits (see table 1). This allows to make confident 
predictions for the probability of observing the value of 
II = 23% which corresponds to the observations of Szalay 
et al. (1993). From (Eh one easily concludes that the value 
of p(23%) ranges in the interval 10~ 8 — 10 -10 . This in- 
terval is determined by taking into account the statistical 
uncertainties of the fits. This implies that the probability 
of obtaining a peak with the observational power-content 
is at most 10~ 8 . Needless to say that this result completely 
discards the possibility of detecting the observed period- 
icity in a universe with this kind of structure. 

As we have stressed, we have performed a similar anal- 
ysis replacing the step functions by Gaussian functions 
with the result that the upper bound for the probability 
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Fig. 5. Power law probability (not normalized) of obtain- 
ing a peak with a corresponding value II. Note that the 
higher the relative power content of the representative 
peak the lower the number of directions of sight where 
such a peaks is present. In particular the probability of 
observing a peak containing 23% of the total integrated 
power is ~ 1CP 8 . 



length scale and no other large scale order should be con- 
sidered as constituting the most viable explanations for 
the observations. 

The conclusion and the main result of our analysis is 
that the existence of a specific scale in the distribution of 
mass density of the universe, as showed by a two-point 
correlation function, does not by itself provide a satisfac- 
tory explanation for the spatial periodicity of 128/i _1 Mpc 
in the galaxy distribution of deep pencil-beam surveys as 
it results in a negligible probability for such observation to 
be obtained in a particular direction. This result seems to 
force us again to look for alternative explanations for the 
observations of periodicity in deep pencil-beam surveys. 
On the other hand, it might be interesting to consider 
variations of this type of model (for exemple, changing the 
intersection rule) to identify the features that lead to an 
enhaced probability for obtaining the desired periodicity. 
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p remains essentially unaffected and all the values of p 
remain within the interval determined by the statistical 
uncertainties. 

It should be emphasized that in this work we are con- 
sidering the possibility of detecting periodicity in a partic- 
ular direction under the sole assumption that there exists 
an intrinsic length scale in the galaxy distribution. Our 
results show that this probability is completely negligible. 
However, this result could drastically be changed within a 
different model which would be characterized not only by 
a length scale but also by the presence of a specific cor- 
relation in the spatial distribution of galaxies. A concrete 
example of a model with this type of additional structure 
has been investigated (Van de Weygaert, 1991; Subba Rao 
& Szalay, 1992). In these works, dynamical Voronoi foams 
(Icke & Van de Weygaert, 1987; Van de Weygaert & Icke, 
1989) with a controlled amount of randomness are pro- 
posed as describing the spatial galactic distribution. The 
probability for finding periodicity in these models depends 
on the amount of randomness and can drastically be in- 
creased from 3% up to 15%. This probability can reach 
even higher values (10% to 40%) in the case of a regu- 
lar lattice (Kurki-Suonio et al., 1990). In the extreme and 
naive model of galaxies distributed on the surfaces of con- 
centric shells, one obviously would obtain a 100% proba- 
bility for an observer at the center. We want to stress that 
the models that are more easily reconciled with acceptable 
cosmological scenarios are those that introduce a priori 
the least amount of scale order, and such that the order 
introduced corresponds to the smallest possible scale. In 
this sense the models that introduce just a characteristic 
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